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The grand canonical ensemble has been used to make predictions for composite yields using 
simple models for nuclear fragmentation. While this gives correct model prediction for high energy 
collisions, it can give very inaccurate results at intermediate energy. 
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A very simple but very popular model for nuclear multifragmentation is this: the nucleus is heated up and breaks up 
into many pieces (composites and new produced particles if the energy is sufficient) strictly according to phase space. 
This occurs in an expanded volume, about three or four times the normal volume. Population strictly according to 
po ' phase-space implies chemical and thermal equilibrium. For simplicity, we will omit new particle production. In such 
cases the number of dissociating particles is fixed. Nonetheless in the past it has been customary for calculational 
simplicity to use the grand canonical ensemble to describe multifragmentation [1-3]. In such a model (GCM), the 
number of particles in the dissociating system is not constant, however one can arrange to have the average number 
_^ , correspond to the actual system. If one is in the classical regime (Fermi or Bose statstics degenerates into Maxwell- 
p,i Boltzmann limit) then the average yield of a composite in the ground state is given by 
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CO , where i is the proton number, j is the neutron number of the composite, fiz is the proton chemical potential, /x„ is 
C~^ ^ the neutron chemical potential and fi j is given by: 
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Here V is the volume within which the particle moves, a = i + j is the mass number of the composite, m is the proton 
mass, g is the spin degeneracy, Eij is the binding energy of the composite and the Maxwell-Boltzmann distribution of 
the momentum of the particle has been integrated over. Usually populations into any states of the composite, ground 
and excited are included (a popular method of including the excited states is to use the Fermi-gas approximation) in 
^ ■ which case fij is replaced by i^ij, the one particle partition function of the particle. Thus we have 



(n,,,)=e^'^^^+^"'3'^"c.,,, (2) 



The Wigner-Seitz approximation of the coulomb energy is usually included [4] and this can be incorporated in the 
?H ' ^i,j by replacing the coulomb self-energy 2:|^ of the composite a(= i + j) to ^^^^(1 — (p/po)^^'^)- The chemical 
potentials /i are fixed from 

^i(71,j) = Z 

where Z, N are the charge and neutron number of the dissociating system. The connection between the GCM and the 
Canonical Model (CM), as described in a textbook, is simple mathematics. Let us denote the CM partition function 
as Q and the GCM partition function as Z, then 



i^hjY 



T-iJ- 



and 



'Corresponding author 



^=n 



y^ ^ii(3ti^-\-jf3fi^)nijQ 



ni, 7=0 



]Jcxp[e(*''^- 



+i/3M,.), 



*j 



10' 



10" 



LJ I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I bd bd I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I bd 



10' 



^10-^ 



>- 



10"^ ^ 



10' 



10" 



T=4.0MeV ^ 




canonical 



^ grandcanonical ^ ^ 



I I I 



T=5.0MeV ^ 




I I I 



20 40 60 80 



20 40 60 80 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I bd ET 



T=7.5 MeV I 




I I I I 111 I I IM 



I ' ' ' I ' ' ' I ' ' ' I 



T=20.0 MeV 




J I I I I I I I I I L 



I . x\ , I 



20 40 60 80 



8 12 16 20 

z 



FIG. 1. The yields in the multifragmentation model using the Grand Canonical ensemble and the canonical ensemble for 
A = 200, Z — 80. Note that the two ensembles give very different values at T=4 MeV. 

Recently [5,6] it has become possible to use the canonical model to calculate yields of fragmentation, whereas, in 
the past, the GCM was universally used. In the general case many composites are present and hence 
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The sum is over all partitions of Z, N into clusters and nucleons subject to two constraints: ^ 
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Si 7 J'^j J — ^- These constraints would appear to make the computation of Qz,n prohibitively difficult which used 
to be the primary reason for the use of the grand canonical ensemble where with two constants ^z , Mn one merely 



arranged the average values to be Z and N . It has been recently realised that a recursion relation exists which allows 
the computation of Qz.n quite easy on the computer even for large Z or N [7]. Three equivalent recursion relations 
exist, any one of which could be used. For example, one such relation is 
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The average number of particles of the species i,j is given by 
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All nuclear properties are contained in Wi 
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FIG. 2. The same as in Fig.l except for the system A = 50, Z = 25. Again note the discrepancies at T=4-0 MeV. 



We are in a position to check, in the nuclear case, the grand canonical predictions for yields with the canonical 
values where particle number conservation is strictly enforced. We show this in Fig.l for A = 200, Z = 80, N = 120 
(on the nuclear scale, a large system) and in Fig. 2 for A = 50, Z = 25, N=25. As there are too many composites. 



we compare isotope yields (yields of the same Z are added up and then compared). The GCM and CM predictions 
are quite close for high temperatures (T >20 MeV) but at low temperatures (w 4 MeV) in spite of 200 being a large 
number, the GCM predictions are significantly different. One encounters such temperatures in intermediate energy 
heavy ion collisions thus one would conclude that one should not use the GCM in intermediate energies. GCM has 
been used at 100 McV/nucleon beam energy in the lab [8] where it may be just adequate but it has also been used at 
much lower energy [9,10] where its usage is very questionable. 

One has to ask if, when the two predictions differ, are both of them wrong or only one of them? If we consider 
thermal and chemical equilibrium to be the fundamental ingredient of the model, then the only way the present 
calculation in CM can be wrong is, if the approximation to quantum statistics (eq.3) is invalid around 5 MeV. 
Following a recent paper [11] we show this is not true and thus CM model results continue to be good. Here is a gist 
of the argument from that paper. 

We use large volumes (3 or 4 times the normal volume). At low temperatures (« 4 MeV) where one might imagine 
the approximation to fail, it survives because different composites appear thus there is not enough of any particular 
species to make (anti)symmetrisation an important issue. At much higher temperature, the number of protons and 
neutrons increase but as is well-known, the n! correction takes the approximate partition function towards the proper 
one. In a different world, the problem could get very difficult. Such a scenario would arise if the physics was such 
that at low temperatures we only had neutrons and protons and no composites. An even worse situation would arise 
if we had only neutrons (or protons). With these preliminaries, let us proceed to to estimate quantitatively the errors 
involved in actual cases that one might encounter. For simplicity only, in this section we will not put in excited states 
of composites and we have not put in the Wigner-Seitz correction although that could have been retained without 
much extra work. 

A recursive relation similar to eq.4 exists even with incorporation of quantum statistics but w^j's are no longer 
one-particle partition functions. We illustrate this first with the example of only protons filling up orbitals i, j, fc... in 
a box. Now 
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The coefficient of e^^^^ is x^ = - — -j, — X^i ^ ^P'^i ^ The canonical partition function for Z protons is very similar to 
eq.4: 
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Here Qo is 1- When the expressions for Xfc are used in the above equation, orbitals are given occupancies greater than 
one and then eliminated by subtraction. This can lead to severe round-off errors when applied to degenerate Fermi 
systems but will not affect the applications here. The number of protons Z is given by 

Z ^ xi—- ^2x2— h Zxz-pr- (8) 

Wz wz ^z 

For generalisation, we will call Xk above as y[ q. The symbol 1,0 means it is a composite with 1 proton and neutron. 

The symbol k means it is obtained from the k — th term in the expansion; y\ q will contribute to Xk.o- 
If instead we had a boson, a deuteron, for example, we would have 

ln[Qgr.can{P, Mp, Mn)] = ^ "^^(^ " e^^.+P^^ ^-P^i) (9) 
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Thus in the case of deuterons y\ \_ (which would contribute to Xk^k) is given by ^^ ^e" "''^ . 

We can treat an assembly of protons, neutrons, deuterons, tritons...etc. The recursive relation if the dissociating 
system has Z protons and N neutrons is 



QzM = y 22 iXijQz-i,N-j (11) 

i=l.ZJ=0,N 

The average number of a composite with ii protons and «2 neutrons is given by 

< n^^,^2 >= yl['a2Qz-ii,N^i2/QzM + ^yf^^^Q z -21^^-21^1 Q ZM + ■■■ (12) 

Unless one is in an extreme degenerate fermi system, one can evaluate the y factors by replacing sums with integration. 
For example, y^"' = ^^^ E* e""^"' where the sum is replaced by J e~"'^''g{e)de = 2^(2^)^/2. Here V is the 
available volume. We have included the proton spin degeneracy; m is the proton mass. For the deuteron, yi [ = 
i / G~''^'^g{e)de. This is 3 x 2'^/2^(2^^3/2^_^ where Eb is the binding energy of the deuteron. It is clear how to 
compute contributions from other composites. 

We test the accuracy of the yields as calculated throughout the main text by comparing with a calculation where 
the complete theory of symmetrisation and antisymmetrisation is used. Subject only to the approximation that 
summation over discrete states has been replaced by an integration over a density of states, the calculation is exact. 
The results are taken from [11]. We take the dissociating system to have Z=25 and N=25. The lowest temperature 
considered is 3 MeV (one might argue that at lower temperature a model of sequential decay is more appropriate). 
The highest temperature shown is 30 MeV. We take a freeze-out volume in which the composites can move freely 
in three times the volume of a normal nucleus with 50 nucleons. Aside from neutrons and protons we allow the 
possibility of composites. Spins and binding energies for deuteron, triton, ^He and ^He are taken from experiments. 
For higher mass composites the binding energy is taken from empirical mass formulas. For fcrmions, spin 1/2 was 
assumed and for bosons spin was assumed. For each Z we take N = Z — 1, Z, and Z + 1. We present in the Table 
I. average yields of protons, neutrons, tritons, '^He, "^He and the sum of yields of all nuclei with charges greater than 
12. Temperature range of 3 to 6 MeV are of interest to many experiments. We also show the results at 30 MeV. The 
CM approximation for composites is seen to be quite good. 

Granting that below a certain temperature, predictions from a grand canonical model gets to be quite erroneous, 
could one predict when it becomes bad and why? The answer to the first part is simple. Usually, the yield < n^ > 
(or < Ua > where a is the mass number) falls with z but below a certain value of the temperature, the yield, after 
falling initially, begins to rise again. If this happens, one must discard the GCM and do a CM. The rise of yields, 
after reaching a minimum, signifies several interesting features in intermediate energy multifragmentation models. In 
percolation and Lattice Gas model [13,14] this signifies the appearence of a percolating cluster. In thermodynamic 
model, the temperature at which this rise, after the minimum, just disappears signifies a first order phase transition 
(in the infinite system and no coulomb limit) [5]. It suggests that at this temperature a large blob of the system, 
usually identified as a liquid, has just disappeared. It has been shown that at the transition temperature, the specific 
heat at constant volume is very different in GCM and CM although they match admirably at higher temperature 
[15]. The reason for the discrepancy is an unusually large fluctuation in the number of particles in the GCM below 
the transition temperature. We have however no fundamental understanding why such large fluctuations appear in 
the GCM below the transition temperature. 
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TABLE I. Comparision of claculations of average yields and E/A, By exact we mean a calculation with proper symmetry. 
Sum over discrete orbitals in a box has been replaced by integration as is the usual practice. 
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